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Abstract 

We study the excitation of harmonic waves in thin elastic samples by a single dislocation 
in arbitrary motion. We consider both screw and edge dislocations that move perpendicularly 
to the surfaces of the layer. In Fourier space the displacement velocity and dynamic stress 
fields generated by the motion of the dislocations are factored as the product of two terms: 
one depends on the motion of the dislocation only, while the other is independent of it, and 
represents the medium's response. The latter term exhibits poles at frequencies that satisfy 
the dispersion relation of the harmonic modes of the plate. In the case of a screw dislocation 
the modes that are excited are a subfamily of the antisymmetric Rayleigh-Lamb modes. For 
an edge dislocation a subfamily of the symmetric Rayleigh-Lamb modes is excited, as well as 
the lowest lying shear mode. The expression corresponding to a uniformly moving screw is 
worked out in detail; it has singular behavior at velocities coincident with the phase velocities 
of the allowed modes. 



1 Introduction 

Recent experimental results [[[], ^, ||, |], |[ [| in fast crack propagation in brittle materials, such 
as glass and plexiglass, have provided evidence that crack propagation is dynamically unstable, 
and bounds that separate different crack behavior have been identified experimentally in terms 
of the crack tip velocity. For example, when the crack tip velocity surpasses 0.5cr, with cr the 
Rayleigh wave velocity, there are strong sound emissions, oscillations of the crack tip velocity, and 
the surfaces left by the crack become quite rough. 

A scenario for a dynamic instability that could explain in a unified way the experimentally 
observed features of acoustic emissions, crack tip velocity oscillations and surface roughness has 
been recently proposed by one of us 0. This scenario is based on continuum elasticity, and it 
incorporates the fact that a crack in motion can be modeled as a continuous, time dependent, 
distribution of dislocations. These moving dislocations, when propagating in a medium that is 
itself oscillating, are subjet to two forces: one is the usual Peach-Koehler force, and the second a 
force that does no work [^J , in this sense analogous to the Lorentz magnetic force acting on a moving 
charged particle. This latter force that does no work is proportional to the dislocation's velocity, 
perpendicular to it, and it is also proportional to the time derivative of the particle displacement 
generated by time dependent external loads. The dynamic instability would occur by the following 
mechanism: 1) the crack begins to propagate and starts emitting elastic radiation; 2) normal modes 
of vibration of the plate are established, and they are responsible for a new oscillating force that 
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does no work and that is proportional to the crack tip velocity; 3) the crack responds to this new 
periodic forcing by emitting radiation mainly at the same frequencies at which it is being forced. 
This radiation reinforces the normal modes, that in turn act on the crack tip: this is the instability. 

The samples of brittle material that are used in experiments jj], ||, || |[ |[ || are typically 
a few mm thick, and ~ 20 cm in width and length. Thus, it appears that useful insights into 
the underlying physics may be obtained by considering the samples as thin layers; that is, elastic 
layers of finite thickness and infinite in-plane dimensions. In this case the elastic sample is a wave 
guide whose possible horizontally polarized shear modes, and its mixed character Rayleigh-Lamb 
modes are well known |^ . In this work we have investigated the excitation of harmonic waves in 
such a wave guide by a dislocation, perpendicular to the surfaces of the layer, in arbitrary motion. 
This represents a first step towards a firm analytic foundation for the scenario proposed in the 
previous paragraph. The principal results of this work are expressions, in Fourier space, that 
describe the displacement velocity field and the stresses due to screw or edge dislocations moving 
perpendicularly to the surfaces of the plate with arbitrary velocity, and the identification of the 
specific wave modes excited in each case. The relevant expressions are given as convolutions of 
a source function localized at the dislocation position with the medium's response, of which the 
Rayleigh-Lamb modes are poles. 

In order to proceed, we use a formulation reported elsewhere 1 1 1 for the stress and velocity 
generated by a moving dislocation in finite elastic samples. This formulation is a generalization 
of line integral representations first found by Mura JlO[ in the case of an infinite medium. The 
specialization of these results for a thin plate is given in Section 2. Section 3 contains the deter- 
mination of amplitudes, whose precise definition is given in the text, associated with each possible 
harmonic mode. Expressions for particle velocity in terms of those amplitudes and a source func- 
tion associated with the moving dislocation are given in Section 4 for a screw and in Section 5 
for an edge dislocation. Section 6 discusses the corresponding stresses. Section 7 carries these 
results further in the special case of a uniformly moving dislocation, with the well known result of 
a dislocation in uniform motion in an infinite medium carried along for comparison. A discussion 
of the implications of these results for our current understanding of the brittle fracture of thin 
plates is offered in Section 8, and Section 9 has concluding remarks. The Fourier components of 
the Green's function of a thin plate |12| are written down explicitely in an Appendix. 



2 Displacement velocity due to a dislocation loop moving 
in a finite sample 

We consider small displacements U m (x, t) in a homogeneous elastic medium of density p and elastic 
constants Cy/y ||. A general formula for the velocity and stress fields generated by a dislocation 
loop undergoing arbitrary motion in an infinite medium was obtained some years ago by Mura JTo| . 
It is written as a convolution of the medium's impulse response with a source localized along the 
dislocation loop (i.e. independent of the slip plane). These formulae have been generalized to the 
case of a bounded elastic medium, as well as the formula for the displacement field that involves 
integration over succesive slip planes The latter expression for the displacement field is: 

U m (x, t) = -bi f dt' f dS' 3 a% n (x', x- t-t'), (1) 
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(attention should be paid to the order of the variables) where a^ m (x, x';t — t') is the elastic stress 
associated with the Green's function: 

a% m {x,x';t-t') ee C ijkl —G lm (x,x';t-t') . (2) 

OXk 

Eq. (Q) is the analog of Mura's relation [lCj ]. Since the problem is homogeneous under time 
translations it is possible to obtain an expression for particle velocity, the time derivative of particle 
displacement U m that involves an integral along the dislocation loop only. Indeed, since dG m k/dt = 
—dGmk/dt', the displacement velocity becomes: 

8TI C°° 8 ( ( 

m -(x,t)=bi dt'{ — [ dS' 3 af™{x\x-t-t')\- dS^ m (x',x;t-t')}. (3) 

J-oo 01 JS(t') J^r(t') 



dt 



The first term is zero since the Green's function vanishes at t—t' = ±oo. Also, since f^sft'j/dt' ^'j = 
Cj pq f L ( t ,)dl' q V p (x' ,t'), where L{t') is the dislocation loop bounding the slip plane S(t'), V p (x\t') 
is the loop's local velocity and tj pq the completely antisymmetric tensor in three dimensions, we 
have: 

^(x,t) = -b i dt' dl' g a? m (x',x;t-t')e jpq V p (x',t'). (4) 

Ci J-oo JL(t') 

This expression for the particle velocity generated by a single dislocation loop within a finite elastic 
sample can also be obtained as a special case of the formulation of Chishko |[2| for the particle 
velocity generated by a continuous distribution of dislocations in a thin elastic plate. 

We now specialize to the case of an infinite plate of thickness 2h, with plane boundaries located 
at z = ±h. The appropriate Green's function was calculated by Chishko jl2) and it has the form: 

G im {z, z'; R-R'-t- t') = GT m (x -x';t- t') - H im (z, z'; R-R';t- t') , (5) 

with the Green's function appropriate for an infinite isotropic medium, and Hi m a kernel that 
cancels the normal stresses due to G°^ at the free surfaces. Hi m is the following convolution: 

(t-t') 



H im (z,z';R- R'-t-t') = Q(t-t') V / dr dS, 

„-i i-JO J 



a=l,2 ' 

xGf„(z, e\ R - R"; t-t'- T)a% m (e - z'- R" - R'; r) , 

(6) 

with d5" = ±dxdy for a = 1, 2 or z" = ±h; 6 is the step function, and £ a = ±h. Here cr^ m is the 
elastic stress associated with the infinite medium Green's function, and Gf n {z, £ a ; R — R';t — t') is 
the Green's function of the homogeneous boundary value problem of the plate: it corresponds to 
the displacement in direction i produced by a point force applied at the surface in direction n at 
position (£ a , R'), and time t': 

P^Gf n {z, T; R-R';t- t') - C iilm -^^G p mn {z, C; R — R';t — t') = , (7) 
with £ a = ±h, and subject to the following boundary conditions at the surfaces: 



_d 

c 

±S ab S m S(R - R')S(t - t') , (8) 



°C n (£ >£°; R — R''it - 1') — c iz i m -g-^G^ in (£, b ,li a ',R-R'',t~t') 
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where ± corresponds to a = 1, 2 or to the surfaces z' = £ a = ±h respectively. 

For convenience, we have reproduced the explicit expressions for the Fourier transforms of 
G? n (z,£;R- R';t-t'), g%£(z,£\k), from Ref. |j in an Appendix: 

gfiZiz, d\k) = f°° dtf dRe-^ t+i ^G^ n (z, £; R; t) . (9) 

3 Harmonic modes of a thin plate 

Harmonic modes in a thin elastic layer have properties, such as polarization and dispersion re- 
lations, that are well known ||. In particular, the dispersion relation can be read off from the 
poles of the Fourier transform of the Green's function Gfj, gf"{z\k) (see Eqs. (|57|)-(|60|) in the 
Appendix). Poles of gf- u '(z\k) occur when A a (k 2 ,ui 2 ) = 0, A s (k 2 ,uj 2 ) = and when sinh(g t /i) = 0, 

cosh(qth) = 0, with qij = y k 2 — uj 2 /cf t , and ci (resp. ct) are the speed of longitudinal (resp. 
transverse) waves in bulk material. 

The poles A a (k 2 , u> 2 ) = of g^ u (z, £|fc) give rise to antisymmetric modes, and we now determine 
their amplitudes, which will be needed in the next section, by taking the limit: 

(t£)i= lim A a9 r»(z,t\k). (10) 

In this way three linearly independent modes, labelled by "j" are obtained. They are vectors 
whose components are labelled by "£" . These antisymmetric modes, of mixed character, are for 
j = f3 = 1, 2 = x, y and j = 3 = z respectively: 

v^(z\k) = ikf)[{k 2 + q 2 ) smh(qih) cosh(g t z) - 2qiq t sinh(g t /i) cosh(qiz)]z 
k ka 

+{-^Tlt[{k 2 + q 2 ) smh.{qih) sinh(g t z) - 2k 2 smh{q t h) smh( qi z)} 



-(Sap TT')[ T J — s'mh(qih) sinh(g t z) 



k a kp (k 2 +q 2 ) 2 
k 2 ' [ q t 
—Aqik 2 ta,nh(q t h) cosh(qih) sinh(q t z)]}e a 

v A (z\k) = qi[{k 2 + q 2 ) cosh(q t h) cosh(q/z) — 2k 2 cosh(qih) cosh(q t z)]z 

-ik a [(k 2 + q 2 ) cosh(q t h) sinh(g/z) - 2q t q t cosh(qih) smh(q t z)]e a 



(11) 



A mode that is not linearly independent of the latter can be obtained as: v A (z\k) = —kpv A (z\k). 
It will be seen below that this is the mode excited by a moving screw dislocation: 

v A (z\k) — —ik 2 [(k 2 + q 2 ) smh(qih) cosh(g t z) — 2qiq t smh(q t h) cosh(qiz)]z 

+ 1tk a e a [2k 2 siah(q t h) sinh(g/z) — (k 2 + q t ) sinh(g;/i) smh(q t z)\ . (12) 

In an analogous way, the amplitude of the symmetric modes corresponding to A s (k 2 ,u> 2 ) = 
are obtained as: 

(4),~ Urn A„gR u (z,t\k). (13) 
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It will also be seen below that the mode excited by a moving edge dislocation corresponds to 
.7 = 3 = z: 

v z s (z\k) = qi[{k 2 + q 2 ) siuh(q t h) smh(qiz) — 2k 2 sinh((fr/i) sinh(q t z)]z 

- ik a [(k 2 + q 2 ) svah(q t h) cosh(g ; z) - 2qiq t sinh(g//i) cosh(g t z)]e a . (14) 

The shear normal modes are of transverse character and polarized parallel to the surfaces of the 
plate. The symmetric shear modes can be obtained by multiplying gf^ by sinh(g t /i) and taking the 
limit smh.{q t h) = 0, or q t h — inn, n — 1, 2, . . .. Thus: 

Vs(z\k) = k 2 (8 aP - ^p^) cos(^z)e Q . (15) 

The antisymmetric shear modes correspond to multiplying gfj(z, £|/c) by cosh(q t h), and then taking 
the limit cosh(q t h) = 0, or q t h — i(2n + 1)tt/2, n = 1, 2, . . .: 

v^(z\k) = k 2 (S af 3 - ^pp) sin(^-^-—^z)e Q . (16) 

Note that the modes with dispersion relation lu = ci.tk that occur in an infinite isotropic medium 
are not normal modes of the plate. The lowest shear mode of the plate has dispersion relation 
u> = c t k, but it is only polarized in the plane of the plate. 



4 Displacement velocity due to a screw dislocation moving 
perpendicularly to the surfaces of the plate 

We consider a straight screw dislocation that is parallel to the z axis, or perpendicular to the 
surfaces of the plate, with a Burgers vector b = bz, in arbitrary motion. The dislocation as a curve 
is described by the location of its points as a function of time: X((t,t) = (Xi(t),X2(t),<t), with 
t the time, and a parametrizing the curve. However, the dislocation loop or integration curve 
L(t') of Eq. should also include segments on the surfaces of the plate and one at infinity (that 
it is easily seen not to contribute), in order to close the loop. The segments on the surfaces of 
the plate do not contribute to the displacement velocity due to the boundary condition satisfied 
by the Green's function on these surfaces: the indices p and q in Eq. (|]) are different from z 
on the segments on the surface, and so, because of the term ej pq , j must be equal to z, thus 
&ji m {x' , x; t — t') is evaluated at the surfaces (x' — x' s ) where it vanishes. The expression in Eq. 
for the displacement velocity is then: 

/oo QX 

dr J daC za kie a p-Q^-Q^jG m k(z, z'\ R — R';t — t) , (17) 

with the convention that greek indices take values (1, 2) or (x, y), and 612 = —£21 = 1, en = £22 = 0. 
Thus: 

dT-g— J da{--^G mz (x 1 X(a,T)- 1 t - t) 
+-^G ma (x,X(a,T);t-T)}. (18) 
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The in-plane Fourier components of the displacement velocity, are defined as: 

y£(z\k)= f dt f dRe~ luJt+llfi V m (R,z,t) . (19) 

J — OO J 

Replacing this into ( |l8|) leads to 

V%(z\k) = -e a ^(k)F^ a (z\k) , (20) 



/■OO _ _ av 

1^%) = / dTe~ WT+lk - x{T) --^{T) , (21) 



with 

/*°° -. _ 

K(h\ = I J TP -^T+Ik-X(T)[ 

dr 

a factor that depends on the velocity of the dislocation, and F^ a (z\k) a factor that represents the 
medium's response: 

F^ a (z\k) = bpc 2 t {ik a f dag^ z (z, a\k) + f da-^-g^ a (z, a\k)} . (22) 

•J — h J — h 

After some algebra, the following expressions for F£ a (z\k) can be obtained: 



mAi) = g^W ^ 1 ' w *' ' (24) 

These factors F^ a (z\k) have poles when A a (fc 2 ,w 2 ) = 0, which gives the dispersion relation 
of the antisymmetric Rayleigh-Lamb modes of the plate. The terms associated with these poles 
are proportional to the Rayleigh-Lamb mode Vj^{z\k) introduced in Eq. (|l2|). This shows that 
the principal contribution (it comes from A a — > 0, or oj ~ the dispersion relation of the 

antisymmetric Rayleigh-Lamb modes) to the displacement velocity field due to a screw dislocation 
moving perpendicularly to the surfaces of the plate comes from the excitation of the antisymmetric 
Rayleigh-Lamb modes associated with (v^) u (z\k) . When Fourier transforming back to configura- 
tion space it becomes clear that only the poles contribute, as they should. An explicitely worked 
out example is given in Section 7. 



5 Displacement velocity due to an edge dislocation moving 
perpendicularly to the surfaces of the plate 

We consider a straight edge dislocation parallel to the z axis, or perpendicular to the surfaces of 
the plate, with a Burgers vector b = b x x + b y y = 6 7 e 7 , and moving arbitrarily. According to Eq. 
(0), for this edge dislocation the displacement velocity becomes: 



dr dz' , G m k(z,z';R — Bf;t-T) 

-oo J-h ° T ax l 



(25) 

R'-X(t) 



Here again, there is no contribution from the segments of the loop L(t') that are on the surfaces of 
the plate and at infinity, for the same reasons that were given for a screw dislocation. Calculating 
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the Fourier components of the displacement velocity, one obtains them as the product of two 
factors: 

V%{z\k) = -e al3 I%{k)PZ a {z\k), (26) 

where Zg(fc) is the factor dependent on the dislocation's velocity, Eq. (pi]), and P^ a (z\k) is a 
factor that represents the medium response: 



It L^s\ K ) 



with: 



A.Ab.L:^ 2 ) = -i^-il[( 7 2 - 2)k 2 b a + 2k a k p b p ] , (29) 



B ia (b\k\u 2 ) = ^ { (7 2 2 ) ba + 2 f| fc b (30) 



" - 2 C t 

C^ OL {b : k,uj ) = i{k a b^ -\- S a ^kpbp 2 — ^ kpbpk&kyj , (31) 

and 7 2 = (q/c 4 ) 2 . Note that the terms of P^ a (z|A;) which have poles at A s (k 2 ,uj 2 ) — are 
proportional to the symmetric normal mode of the plate, (vg) u (z\k) , of Eq. (|l4|). There is also a 
pole at g 2 = with non zero residue associated. This shows that the principal contribution (it comes 
from A s — » 0, or to ~ (fc): the dispersion relation of the symmetric Rayleigh-Lamb modes) to 
the displacement velocity comes from the excitation of the Rayleigh-Lamb mode (vg)^ '(z\k) , and 
also from the lowest lying shear mode. As in the case of the screw dislocation, going back to 
configuration space shows, via Cauchy's theorem, that only those modes contribute. 

6 Dynamic stresses on the plate 

Since the elastic stresses satisfy: 

Q Qv m 

g^Cij = Ciji m -Q^-(x,t) , (32) 



for lo 7^ its Fourier components are: 

a%(z\k) = ^V l V£(z\k), (33) 



with Vf = {—ik x , —iky, d/dz). We shall check that this expression satisfies the appropriate bound- 
ary conditions at the free surfaces. 

6.1 Screw dislocation 



Using Eqs. (p0|), (|33|), one obtains the following expression for the stress produced by a dislocation 
that moves perpendicularly to the surfaces of the plate (u ^ 0): 

a^ z (z\k) = ~I^(k)e a0 K^ z (z\k) , (34) 
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with: 

KZ(z\k) ee pcU^k 1 { 1 2 ~2)F- a (z\k)+l 2 ^ a {z\k)} 

k%{z\%) ee pc 2 {-ik s F? a (z\k) + -^F? a (z\k)} . (35) 

Using the expressions obtained for F£ a (z\k) (see Eqs. (p3|), (E4J)), one obtains: 



K%?(z\k) ee bpc 2 ^{\ - —[(k 2 + q 2 ) 2 smh(q l h)cosh{qtz)-4qiq t k 2 smh(q t h)cosh{q l z)}} . (36) 



>kgk a 1 

If one evaluates the stress on the surfaces of the plate z = ±h, one sees, using Eq. (^), that 
it vanishes, as it should for a traction free surface. Also, 

K™(z\k) = bpe$ 2thak ^ +9 ^ {sinh( gt /t) sinhfez) - sinh(g^) sinhfez)} . (37) 

This is also zero on the surfaces of the plate z = ±h, as it should. The other components of the 
stress also decompose as the product of two factors. 



6.2 Edge dislocation 

The Fourier components of the stress, a^ z (z\k), produced by an edge dislocation moving perpen- 
dicularly to the surfaces of the plate (lj ^ 0) are decomposed also into two factors: 

<Mk) = ~e a ^(k)Q^ z (z\k) , (38) 



with: 



and: 



QTAAk) = D a (k, b){l - — 2 [(k 2 + <fi) 2 smh(q t h)cosh( qi z) 



4k 2 qiq t sinh(gj/i) cosh(g t 2:)]} , (39) 

(40) 



_ 2ik (fc 2 + q 2 ) -> -> 

Q"%(z\k) = — ~klTa — ir D a(k, b){smh(q t h) sinh^/z) - smh(qih) sinh(g t z)} , 



D a (k, b) ee pc 2 ^ 2) [( 7 2 - 2)k 2 b a + 2k a k p b p ] . (41) 

7 



These expressions for Q™J(z|fc) assure that the components of the stress af z are null on the surfaces 
of the plate. 



7 Example: a uniformly moving screw dislocation 

We consider the case of a screw dislocation that moves perpendicularly to the surfaces of a plate, 
with uniform velocity V a x. We will make a parallel with the case of the same dislocation moving 
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in an infinite medium, for which the particle velocity is well known |13|| . In the case of a thin layer, 
according to Eqs. (|l9|), (|20|), the displacement velocity in coordinate space is: 



V m (x,t) = -J^^p JdkJ du;e^ k - s I^k)F^ a (z\k) . (42) 

For a screw dislocation moving in a plate the factor F^ a (z\k) was given in Eqs. (|23|), (|2~i|). In an 
infinite medium the analogous quantity is: 

= , (43) 

which shows that, in this case, only shear modes polarized in the z direction are excited, since 
there is a pole at g t 2 = or w = ±c t k. Notice that in the case of the plate there is an "apparent" 
pole at qf = also, but its residue is zero, meaning that "infinite medium" shear modes are not 
excited. 

Using the definition of Ip(k) given in Eq. (|2l]), Eq. ( p2[ ) for the displacement velocity field 
component m = 3 = z in the plate, call it Vf(x, t), can be written as (we will concentrate on the 
component m = 3 = z, which is the only nonvanishing component in the infinite medium case): 

Vffrt) = ~ h ^e a , J dkj to J dTe ^-r)-M S -v^K {l _ . (44) 

The integrand in this expression has poles at A a (k 2 ,oj 2 ) — 0. These poles give the dispersion 
relation of the antisymmetric Rayleigh-Lamb modes of the plate, w = w^(fc) + «£• A small positive 
imaginary part has been added in order to satisfy causality: if r > t there are no contributions 
because one closes the contour of the to integral in the lower half plane, where there are no poles. 
When t > t we do the integral over uj by closing the contour on the upper half plane, and one 
obtains the contribution from all the poles already mentioned: 



bV _ d ^ f j£ f j iu n {k)(t-r)-ik-(x-V Tx) _ ( v a)z i z \k) 



(27r) 2 «* dXa Z^J q 2d -t{k\^) 



, (45) 

LO — UJ n (k) 



while the analogous expression in the infinite medium case for the z component of the displacement 
velocity field, V£°(x, t), comes from calculating the residues at the poles w n = ±ctk. This gives: 

ihV r 2 e: fl r _ /"* - p ikct(t—r) p -ikct(t-r) 

vrm - -^tJ s L j***- y " "(Ssr - ssr* ■ (46) 

The last two equations show that the solution for V z (x, t) can be written as a superposition of 
contributions from different modes. The integration over r can be done, and it renders for the case 
of the plate: 



ibV d ^ f t e -Mz-v t±) (v%)<2(z\k) 



(2tt)2 dy^J K(fc) - k x V Q ) g t 2 £A a (fc2,^) 

and in the case of the infinite medium: 

bV c t d f dk ik.^-Vatx), 1 1 



(47) 



.(fe) 



V ™ {x ' t] 2(27r) 2 9yi fc 6 * ° " ^ (fcc t - k x V ) {-h* - k x V Q ) ] (48) 
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These last two equations show that there would be a large contribution to the final result from the 
region where the "Doppler shifted" frequency of the modes {ui n {k) — k x V Q ) is approximately zero 
(if that is possible for a given V ). For example, in the infinite medium case if the velocity V of the 
dislocation goes beyond Ct, one gets Cerenkov type radiation. But, for the case of the plate, the 
phase velocities u) n (k)/k x are always greater than Vr, the Rayleigh wave velocity, meaning that 
there is no large excitation of modes at uniform velocities V a < Vr- 

The angular integral over k can't be done analytically for the case of the plate. We continue with 
the infinite medium case in order to illustrate that one can get the exact solution as a superposition 
over the shear modes excited. If k — k(cos8, sin (9), then Eq. ( fl8|) can be written: 

ibr Z" 00 r* 11 
V ^ ^| * /^c»B i a^--^«^*» { ___ + j^^} . (49) 

Doing the integral over k: 

bc t f n ,„ sin8 1 1 

z ^ X,t > = 2(2^)2 y_ w ((x-V o t)cos6 + ysm0y(%-cos0) + (f + cos0) } ' ( } 

This angular integral can be done by using the change of variables z — e t9 , then it becomes an 
integral over the unit circle in the complex plane, and then one applies Cauchy's theorem. The 
result is: 



v.m-^ ^'' l " )2 , . (5i) 

which corresponds to the exact solution that can be obtained in a simpler way |l3[ : 



b ,Wi-(f) 2 

^t) = -tan-H ), (52) 

This is, of course, the static solution for a screw dislocation in an infinite medium, U z = (b/2ir)(f> = 
(b/2ir)tan~ 1 (y/x), when Lorentz-transformed to a frame moving with velocity Vq. 



8 Relation with experiments on fast fracture 

Experiments on fast fracture of glass and plexiglass Q, ||, Q have been made in mode I 
loading, which means that, in order to use our results, the crack should be modeled by edge dislo- 
cations. Our results about the motion of edge dislocations show that the mode excited in this case 
is the Rayleigh-Lamb mode associated with (vg) u (z\k). Moreover, experiments exhibit acoustic 
and surface roughness spectra that have strong peaks at frequencies corresponding to the lowest 
lying long wavelength Rayleigh-Lamb modes |7| , and it is tempting to think that long wavelength 
modes, requiring less energy, will be more efficiently excited. The long wavelength solution for the 
dispersion relation of these modes comes from taking the limit k — > in A s (k 2 , ui 2 ) = 0, and leads 
to two possibilities: cos(ujh/ci) — or "odd" frequencies u n = (2n + l)irci/2h (n = 1,2, .. .) and 
sm(ujh/ct) — or "even" frequencies u> m — mirct/h (m = 1,2,.. .). The limit k — » for the mode 
v z s (z\k) is: 

3 

v%{z\k) ~ — i — ^ sin( — h) sin( — z)z (53) 

ClCf C t Cl 
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This expression goes to zero for the frequencies previously called "even" , which means that the 
amplitude of excitation of these modes will be lower than the corresponding amplitudes associated 
with the "odd" frequencies. The fact that experiments deal with a crack, that is with many 
dislocations side by side, rather than a single dislocation will not change this fact and it may well 
happen that deviations from straight-front crack tip behaviour must be considered. 

9 Concluding remarks 

In this paper we have given an expression for the elastodynamic fields that arc excited in a thin 
plate by a dislocation, perpendicular to the faces of the plate, in arbitray motion. This expression 
is a convolution of the impulse response of the plate with a source function that is localized at 
the dislocation position. An explicit expression for the impulse response has been determined in 
Fourier space as a rational function, with poles at the well known Raylcigh Lamb modes, both for 
screw and edge dislocations. Application to the case of a uniformly moving screw shows that large 
amounts of radiation are to be expected when the dislocation's velocity coincides with the phase 
velocity of one of the modes. 

This work is part of an ongoing project that attempts to understand recently observed dynamic 
crack instabilities in thin plates in terms of the interaction of the crack tip with the wave modes 
of the plate. We hope the formalism presented in this work will prove helpful in advancing this 
project. 
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Appendix: Fourier components of the plate's Green's func- 
tion 

In this Appendix we reproduce the expression for the Green's function of a thin plate given by 
Chishko |L2|] . According to Eq. (||) the Fourier components of the plate's Green's function are: 

gf m (z, z'\k) = g°^{z - z'\k) - h^(z, z'\k) . (54) 

From Eq. (JgJ) : 

w m (z,z'\k) = ^(±)[-%c„ z ^ w (^,ri^)^(r-^|fc) 

a 

+ C nzzp g?»(z,e\k)^g™(e-z'\k)], (55) 

where (±) refers to a = 1, 2 or £ a = ±/i. The Fourier transforms of the infinite medium's Green's 
function are Jl2| : 

1 , , p-1i\ z - z '\ in 2 p-1t\z-z'\ 

g°°"(z z'\k) = --^{s 4 Wj + 72 } - (56) 

with s*'* = {ik x ,ik y ,qi t ts{z — z')), s{x) = sign(x), and q^t = k 2 — uj 2 /cf t . The Fourier transforms 
of the Green's function (see Eqs. (@), (|)) are ||: 

9^(z,Z\k) = AE M^{z,d\k) . (57) 
pCt \=i.t 



Here M^(z,£\k) are matrices with components: 

M l z %= qi (k 2 + q?)$ 2 (z\0 ■ 

M% = -2iqiq t k fj ^ A {z\0 

M l £ = ~l{k 2 + q 2 t )k a <S>M0 
M l ^ = -2q t k a k ^ x {z\£) 



M% = -2 m k 2 ^\z) 
M% = l(k 2 + q 2 )k^\z) 
M* a % = 2iq iqt k a 3> 4 (Z\z) 



12 



™% = ^*(fc 2 + q f)^m + (s a0 - ^)[ (fc + qt Qt) ±mk 2 m\z)\ 

where £ is equal to either h or — h. The functions $ p (u\v) and &(u\v) are defined as: 
$>i(u\v) = -r— cosh(g;u) cosh(g t w) + — — sinh(gju) sinh(<7 t u) 
$>2{u\v) = sinh(<7;M) sinh(g t w) + cosh(<7/u) cosh(q t w) 

$s(u\v) = — — cosh(g;u) sinh(<7i«) + — — sinh(<ftti) cosh(q t u) 
$>4(u\v) = sinh(<7;ii) cosh(qiw) + cosh(gju) sinh(<ftu) 

fy(u\v) = coth(q t u) sinh(qiu) cosh(q t v) 

+ —— tanh(gtu) cosh(g;u) smh(qtv) , 

with: 

A s (k 2 , uj 2 ) = (k 2 + q 2 ) 2 cosh(qih) smh(q t h) — 4k 2 qiq t cosh(q t h) sinh(qih) 
A a (k 2 ,uj 2 ) = (k 2 + q 2 ) 2 smh(qih) cosh(q t h) - 4k 2 qiq t smh(q t h) cosh(<?;/i) 
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